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Abstract. Cosmological consequences of the brane world scenario are reviewed in a 
pedagogical manner. According to the brane world idea, the standard model particles 
are confined on a hyper-surface (a so-called brane), which is embedded in a higher- 
dimensional spacetime (the so-called bulk). We begin our review with the simplest 
consistent brane world model: a single brane embedded in a five-dimensional Anti-de 
Sitter space-time. Then we include a scalar field in the bulk and discuss in detail the 
difference with the Anti-de Sitter case. The geometry of the bulk space-time is also 
analysed in some depth. Finally, we investigate the cosmology of a system with two 
branes and a bulk scalar field. We comment on brane collisions and summarize some 
open problems of brane world cosmology. 
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1. Introduction 

The idea of extra dimensions was proposed in the early twentieth century by Nordstrom 
and a few years later by Kaluza and Klein P . It has reemerged over the years in theories 
combining the principles of quantum mechanics and relativity. In particular theories 
based on supersymmetry, especially superstring theories, are naturally expressed in 
more than four dimensions [2.. Four dimensional physics is retrieved by Kaluza-Klein 
reduction, i.e compactifying on a manifold of small size, typically much smaller than 
the size of an atomic nucleus. 

Recent developments in string theory and its extension M-theory have suggested 
another approach to compactify extra spatial dimensions. According to these 
developments, the standard model particles are confined on a hypersurface (called brane) 
embedded in a higher dimensional space (called bulk). Only gravity and other exotic 
matter such as the dilaton can propagate in the bulk. Our universe may be such a 
brane-like object. This idea was originally motivated phenomenologically (see [H] [I]) 
and later revived in string theory. Within the brane world scenario, constraints on the 
size of extra dimensions become weaker, because the standard model particles propagate 
only in three spatial dimensions. Newton's law of gravity, however, is sensitive to the 
presence of extra-dimensions. Gravity is being tested only on scales larger than a tenth 
of a millimeter and possible deviations below that scale can be envisaged. 

From the string theory point of view, brane worlds of the kind discussed in this 
review spring from a model suggested by Horava and Witten |Hj. The strong coupling 
limit of the E 8 x E 8 heterotic string theory at low energy is described by eleven 
dimensional supergravity with the eleventh dimension compactified on an orbifold with 
Z2 symmetry, i.e. an interval. The two boundaries of spacetime (i.e. the orbifold fixed 
points) are 10-dimensional planes, on which gauge theories (with the E 8 gauge groups) 
are confined. Later Witten argued that 6 of the 11 dimensions can be consistently 
compactified on a Calabi-Yau threefold and that the size of the Calabi-Yau manifold 
can be substantially smaller than the space between the two boundary branes jH] . Thus, 
in that limit space-time looks five-dimensonal with four dimensional boundary branes 
[TU] . This provides the underlying picture for many brane world models proposed so 
far. 

Another important ingredient was put forward by Arkani-Hamed, Dimopoulos and 
Dvali (ADD), [TT] and [T2], following an earlier idea by Antoniadis [THj, who suggested 
that by confining the standard model particle on a brane the extra dimensions can be 
larger than previously anticipated. They considered a flat bulk geometry in (4 + d)- 
dimensions, in which d dimensions are compact with radius R (toroidal topology). The 
four-dimensional Planck mass Mp and the (4 + <i)-dimensional Planck mass Mf un( i, the 
gravitational scale of the extra dimensional theory, are related by 

Ml = m££b*. (l) 

Gravity deviates from Newton's law only on scales smaller than R. Since gravity is 
tested only down to sizes of around a millimeter, R could be as large as a fraction of a 
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millimeter. 



ADD assumed that the bulk geometry is flat. Considerable progress was made 
by Randall and Sundrum, who considered non-flat, i.e. warped bulk geometries [Hj . 
[To] . In their models, the bulk spacetime is a slice of Anti-de Sitter spacetime, i.e. 
a space-time with a negative cosmological constant. Their discovery was that, due to 
the curvature of the bulk space time, Newton's law of gravity can be obtained on the 
brane of positive tension embedded in an infinite extra-dimension. Small corrections to 
Newton's law are generated and constrain the possible scales in the model to be smaller 
than a millimetre. 

They also proposed a two-brane model in which the hierarchy problem, i.e. the 
large discrepancy between the Planck scale at 10 19 GeV and the electroweak scale at 
100 GeV, can be addressed. The large hierarchy is due to the highly curved AdS 
background which implies a large gravitational red-shift between energy scale on the 
two branes. In this scenario, the standard model particles are confined on a brane with 
negative tension sitting at y = r c , whereas a positive tension brane is located at y — 0. 
The large hierarchy is generated by the appropriate inter-brane distance, i.e. the radion. 
It can be shown that the Planck mass Mpi measured on the negative tension brane is 



where M5 is the five-dimensional Planck mass and A5 the (negative) cosmological 
constant in the bulk. Thus, we see that, if M 5 is not very far from the electroweak 
scale My/ ~TeV, we need kr c pa 50, in order to generate a large Planck mass on our 
brane. Hence, by tuning the radius r c of the extra dimension to a reasonable value, one 
can obtain a very large hierarchy between the weak and the Planck scale. Of course, 
a complete realization of this mechanism requires an explanation for such a value of 
the radion. In other words, the radion needs to be stabilized at a certain value. The 
stabilization mechanism is not thoroughly understood, though models with a bulk scalar 
field have been proposed and have the required properties [T5] . 

Another puzzle which might be addressed with brane models is the cosmological 
constant problem. One may invoke an extra dimensional origin for the apparent (almost) 
vanishing of the cosmological constant. The self-tuning idea advocates that the 
energy density on our brane does not lead to a large curvature of our universe. On the 
contrary, the extra dimension becomes highly curved, preserving a flat Minkowski brane 
with apparent vanishing cosmological constant. Unfortunately, the simplest realization 
of this mechanism with a bulk scalar field fails due to the presence of a naked singularity 
in the bulk. This singularity can be shielded by a second brane whose tension has to be 
fine-tuned with the original brane tension ^H|- In a sense, the fine tuning problem of 
the cosmological constant reappears through the extra dimensional back-door. 

Finally, we will later discuss in some detail another spectacular consequence of 
brane cosmology, namely the possible modification to the Friedmann equation at very 
high energy This effect was first recognised in [20] in the context of inflatonary 




(2) 
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solutions. As we will see, Friedmann's equation has, for the Randall-Sundrum model, 
the form (|2U and [22]) 

# 2 = |p 2 + ^P + A, (3) 

relating the expansion rate of the brane H to the (brane) matter density p and the 
(effective) cosmological constant A. The cosmological constant can be tuned to zero by 
an appropriate choice of the brane tension and bulk cosmological constant, as in the 
Randall-Sundrum case. Notice that at high energies, for which 

967rGjv ,.s 
P^— — ( 4 ) 

where k\ is the five dimensional gravitational constant, the Hubble rate becomes 

H oc p, (5) 

while in ordinary cosmology H oc ^fp. The latter case is retrieved at low energy, i.e. 
9QttG n 

P< 1—, (6) 



Of course modifications to the Hubble rate can only be significant before nucleosynthesis. 
They may have drastic consequences on early universe phenomena such as inflation. 

In this article we will review these and other aspects of the brane world idea in 
the context of cosmology. In order to give a pedagogical introduction to the subject, 
we will follow a phenomenological approach and start with the simplest model, i.e. the 
Randall-Sundrum model, with a brane embedded in a five-dimensional vacuum bulk 
spacetime (section 2). Later we will include a bulk scalar field (section 3). In section 4 
we will discuss the geometry of the bulk-spacetime in some detail. In the last part we 
will discuss more realistic models with two branes and bulk scalar fields (section 5). In 
section 6 we will discuss brane collisions. Open questions are summarized in section 7. 
We would like to mention other review articles on brane worlds and cosmology, taking 
different approaches from the one taken here [2Sj-[2Hl- We will mostly be concerned 
with the case, in which the bulk space-time is five-dimensional. 



2. The Randall-Sundrum Brane World 



Originally, Randall and Sundrum suggested a two-brane scenario in five dimensions 
with a highly curved bulk geometry as an explanation for the large hierarchy between 
the Planck scale and the electroweak energy-scale [H|. In this scenario, the standard 
model particles live on a brane with (constant) negative tension, whereas the bulk is a 
slice of Anti-de Sitter (AdS) spacetime , i.e. a space-time with a negative cosmological 
constant. In the bulk there is another brane with positive tension. This is the so- 
called Randall-Sundrum I (RSI) model. Analysing the solution of Einstein's equation 
on the positive tension brane and sending the negative tension brane to infinity, an 
observer confined to the positive tension brane recovers Newton's law if the curvature 
scale of the AdS is smaller than a millimeter [E]. The higher-dimensional space is 
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non-compact, which must be contrasted with the Kaluza-Klein mechanism, where all 
extra-dimensional degrees of freedom are compact. This one-brane model, on which 
we will concentrate in this section, is the so-called Randall-Sundrum II (RSII) model. 
It was shown, there is a continuum of Kaluza-Klein modes for the gravitational field, 
contrasting with the discrete spectrum if the extra dimension is periodic. This leads 
to a correction to the force between two static masses on the brane. To be specific, it 
was shown that the potential energy between two point masses confined on the brane is 
given by 

-r t / s G N mim 2 A I 2 „, /rTX 

v ( r ) = — : — 1 + 3 + °( r ) • ( 7 ) 



r 

In this equation, / is related to the five-dimensional bulk cosmological constant A5 by 
I 2 = —6/(k 2 A^) and is therefore a measure of the curvature scale of the bulk spacetime. 
Gravitational experiments show no deviation from Newton's law of gravity on length 
scales larger than a millimeter [29J. Thus, I has to be smaller than that length scale. 

The static solution of the Randall and Sundrum model can be obtained as follows: 
The total action consists of the Einstein-Hilbert action and the brane action, which in 
the Randall-Sundrum model have the form 

Seh = - / dx 5 ^g~^ ^ + A 5 ) , (8) 

S bra nc= jdx A ^^){~a). (9) 

The parameter A5 (the bulk cosmological constant) and a (the brane tension) are 
constant and k 5 is the five-dimensional gravitational coupling constant. The brane 
is located at y = and we assume a Z 2 symmetry, i.e. we identify y with —y. The 
ansatz for the metric is 

ds 2 = e- 2K ^r)^ u dx tx dx v + dy 2 . (10) 
Einstein's equations, derived from the action above, give two independent equations: 

GK' 2 = - k 2 5 A 5 

3K" = 4ad(y). 
The first equation can be easily solved: 

K = K(y) = \/-yA 5 y = ky, (11) 

which tells us that A 5 must be negative. If we integrate the second equation from — e 
to +e, take the limit e — > and make use of the Z 2 -symmetry, we get 

6K'\ = K 2 5 a (12) 

Together with eq. (fTTj) this tells us that 

A 5 = -^o 2 (13) 
o 

Thus, there must be a fine-tuning between the brane tension and the bulk cosmological 
constant for static solutions to exist. In this section we will discuss the cosmology of 
this model in detail. 
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2.1. Einstein's equations on the brane 

There are two ways of deriving the cosmological equations and we will describe both 
of them below. The first one is rather simple and makes use of the bulk equations 
only. The second method uses the geometrical relationship between four-dimensoinal 
and five-dimensional quantities. We begin with the simpler method. 



2.1.1. Friedmann's equation from five- dimensional Einstein equations In the following 
subsection we will set k§ ee 1. We write the bulk metric as follows: 

ds 2 = a 2 b 2 (dt 2 - dy 2 ) - a 2 5 lJ dx i dx j . (14) 

This metric is consistent with homogeneity and isotropy on the brane located at y — 0. 
The functions a and b are functions of t and y only. Furthermore, we have assumed 
flat spatial sections, it is straightforward to include a spatial curvature. For this metric, 
Einstein equations in the bulk read: 



a 2 b 2 G\ = 
a 2 b 2 G\ = 
a 2 b 2 G\ ee 
a 2 b 2 G\ ee 



a/ 



ab a" 
ab a 

,/2 



ab 



+ kb 2 



2 1,2 

a b 



d ab a' a'V , 9 

h 2 k ^ 

a ao er ao 

do' a'6 



a aa 

— + 2— + — + — 



a 



a" 



ao ao 



= -a 2 o 2 T 5 ° 



Pb + P%-2/&)J (15) 
a 2 6 2 T| (16) 

(17) 



'abb 2 a" b" b' 2 Il2 
a b b z a b b z 



-a z b 



2i,2 



p B +pS(y - y b ) 



3i 



(18) 



where the bulk energy-momentum tensor T fe a has been kept general here. For the 
Randall-Sundrum model we will now take Pb = ~Pb — A 5 and T 5 ° = 0. Later, in 
the next section, we will use these equations to derive Friedmann's equation with a bulk 
scalar field. In the equations above, a dot represents the derivative with respect to t 
and a prime a derivative with respect to y. 

Let us integrate the 00-component over y from — e to e and use the fact that 

a (y) = a (-y)i Kv) = K~y)i a '(y) = ~ a (-y) and b '(y) = -K-v) ( Le - ^-symmetry). 

Then, taking the limit e->0we get 
a' , 1 



\y=Q 



6 



(19) 



(20) 



Integrating the ij-component in the same way and using the last equation gives 

- \y=o = --ab(p + p). 

These two conditions are called the junction conditions. The other components of the 
Einstein equations should be compatible with these conditions. It is not difficult to 
show that the restriction of the 05 component to y = leads to 



P + S-(p + p) =0, 
a 



(21) 
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where we have made use of the junction conditions (|19|) and (|2(J|). This is nothing but 
matter conservation on the brane. 

Proceeding in the same way with the 55-component gives 



a ab , , 9 a 2 b 2 \ 1 
r + kb 2 



—p (p + 3p) + q B 



(22) 



a ab 3 

Changing to cosmic time dr = abdt, writing a = exp(a(t)) and using the energy 
conservation gives ([201, Ell) 

d(H 2 e ia ) _ 
da 

In this equation aH = da/ dr. This equation can easily be integrated to give 

H 2 = ^ + ^ + ^. (24) 




(23) 



The final step is to split the total energy-density and pressure into parts coming from 
matter and brane tension, i.e. to write p = pu + er and p = pu — o~. Then we find 
Friedmann's equation 

A4 £ 

3 a 4 

where we have made the identification 
87rG a 



R 2 - 8nG n 



x PM 

2a 



+ £ + 4> ( 25 ) 



3 18 
A 4 a 2 A 



(26) 



4 O iVc; . _ N 

— = — + — ■ 27 

3 36 6 v ; 

Comparing the last equation with the fine-tuning relation (JT3Jl in the static Randall- 
Sundrum solution, we see that A 4 = in this case. If there is a small mismatch between 
the brane tension and the five-dimensional cosmological constant, then an effective four- 
dimensional cosmological constant is generated. Another important point is that the 
four-dimensional Newton constant is directly related to the brane tension. The constant 
p appears in the derivation above as an integration constant. The term including p is 
called the dark radiation term (see e.g. [S2]-|S11)- The parameter p can be obtained 
from a full analysis of the bulk equations |H2]-[III1 (we will discuss this in section 4). 
An extended version of Birkhoff's theorem tells us that if the bulk spacetime is AdS, 
this constant is zero [HE]- If th e bulk is AdS-Schwarzschild instead, p is non-zero but a 
measure of the mass of the bulk black hole. In the following we will assume that p = 
and A 4 = 0. 

The most important change in Friedmann's equation compared to the usual four- 
dimensional form is the appearance of a term proportional to p 2 . It tells us that if the 
matter energy density is much larger than the brane tension, i.e. pu 3> &, the expansion 
rate H is proportional pm, instead of J~Pm ■ The expansion rate is, in this regime, larger 
in this brane world scenario. Only in the limit where the brane tension is much larger 
than the matter energy density, the usual behaviour H cx Jp~M is recovered. This is the 
most important change in brane world scenarios. It is quite generic and not restricted 
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to the Randall-Sundrum brane world model. From Friedmann's equation and from the 
energy-conservation equation we can derive Raychaudhuri's equation: 



^7 = -4ttG(p m + p M ) 



1 + P_M 

a 



(2? 



We will use these equations later in order to investigate inflation driven by a scalar field 
confined on the brane. 

Notice that at the time of nucleosythesis the brane world corrections in Friedmann's 
equation must be negligible, otherwise the expansion rate is modified and the results 
for the abundances of the light elements are unacceptably changed. This implies that 
a > (lMeV) 4 . Note, however, that a much stronger constraint arises from current 
tests for deviation from Newton's law [39J (assuming the Randall-Sundrum fine-tuning 
relation (fTS]) ): 3 > 10 5 TeV and o > (lOOGeV) 4 . Similarily, cosmology constrains the 
amount of dark radiation. It has been shown that the energy density in dark radiation 
can at most be 10 percent of the energy density in photons 



2.1.2. Another derivation of Einstein's equation There is a more powerful way of 
deriving Einstein's equation on the brane [4*0] . Consider an arbitrary (3+1) dimensional 
hypersurface Ai with unit normal vector n a embedded in a 5 dimensional spacetime. 
The induced metric and the extrinsic curvature of the hypersurface are defined as 

h\ = 8 a b -n a n b , (29) 
K ab = h a c h b d V c n d . (30) 

For the derivation we need three equations, two of them relate four-dimensional 
quantities constructed from h ab to full five-dimensional quantities constructed from g ab . 
We just state these equations here and refer to jH] for the derivation of these equations. 
The first equation is the Gauss equation, which reads 

R ( ah = hJh b k h c l h d m R,Mm - 2K a[c K d]b . (31) 

This equation relates the four-dimensional curvature tensor R a A bcd , constructed from h ab} 
to the five-dimensional one and K ab . The next equation is the Codazzi equation, which 
relates K ab , n a and the five-dimensional Ricci tensor: 

Vi 4) n - Vi 4 >K = rfh\R bc . (32) 

One decomposes the five-dimensional curvature tensor R a bcd into the Weyl-tensor C abcd 
and the Ricci tensor: 

2 1 

Rabcd = g (g a [cRd]b — gb[cRd]a^) — -^Rg a [cgb]d + C a b Cd . (33) 

If one substitutes the last equation into the Gauss equation and constructs the four- 
dimensional Einstein tensor, one obtains 

G { 2 = I (G cd h c a h\ + (G cd n c n d - h^j h ab 

+ KK ab - K a c K bc -\{K 2 - K cd K cd ) h ab - E ab , (34) 
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where 

E ab = C acbd n c n d . (35) 

We would like to emphasize that this equation holds for any hypersurface. If one 
considers a hypersurface with energy momentum tensor T ab , then there exists a 
relationship between K ab and T ab (T is the trace of T ab ) [32] : 

[K ab ] = -4 (r ab - h ab T) , (36) 



3 

where [...] denotes the jump: 

[f](y) = \un^ (f(y + e)-f(y-e)). (37) 

These equations are called junction conditions and are equivalent in the cosmological 
context to the junction conditions ([Hty and ([20]). Splitting T ab = r ab — ah ab and inserting 
the junction condition into equation (|34p. we obtain Einstein's equation on the brane: 

G$ = 87rGr a6 - A 4 /i ab + ^ - £ afe . (38) 

The tensor 7i ab is defined as 

TTab = -T^TTab ~ tV(, C + -h ab T cd T Cd - —T 2 h ab , (39) 

iz 4 o z4 



SttG = (40) 




whereas 



A 4 y ( A,, - ^^ ) . (Jl) 

Note that in the Randall-Sundrum case we have A 4 = due to the fine-tuning between 
the brane tension and the bulk cosmological constant. Moreover E ab = as the Weyl- 
tensor vanishes for an AdS spacetime. In general, the energy conservation and the 
Bianchi identities imply that 

4^ab = V a E ab (42) 

on the brane. 

Clearly, this method is powerful, as it does not assume homogeneity and isotropy 
nor does it assume the bulk to be AdS. In the case of an AdS bulk and a Friedmann- 
Robertson walker brane, the previous equations reduce to the Friedmann equation and 
Raychaudhuri equation derived earlier. However, the set of equations on the brane are 
not closed in general [13], as we will see in the next section. 



2.2. Slow-roll inflation on the brane 

We have seen that the Friedmann equation on a brane is drastically modified at high 
energy where the p 2 terms dominate. As a result the early universe cosmology on branes 
tends to be different from standard 4d cosmology. In that vein it seems natural to look 
for brane effects on early universe phenomena such as inflation (see in particular 
and 45J) and on phase-transitions [4"E] . 
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The energy density and the pressure of a scalar field are given by 



(43) 



(44) 



where V(<f>) is the potential energy of the scalar field. The full evolution of the scalar 
field is described by the (modified) Friedmann equation, the Klein-Gordon equation 
and the Raychaudhuri equation. 

We will assume that the field is in a slow-roll regime, the evolution of the fields is 
governed by (from now on a dot stands for a derivative with respect to cosmic time) 



It is not difficult to show that these equations imply that the slow-roll parameter are 
given by 



The modifications to General Relativity are contained in the square brackets of these 
expressions. They imply that for a given potential and given initial conditions for the 
scalar field the slow-roll parameters are suppressed compared to the predictions made 
in General Relativity. In other words, brane world effects ease slow-roll inflation [44]. 
In the limit a ^ V the parameter are heavily suppressed. It implies that steeper 
potentials can be used to drive slow-roll inflation [32] • Let us discuss the implications 
for cosmological perturbations. 

According to Einstein's equation (|58j) . perturbations in the metric are sourced not 
only by matter perturbations but also by perturbations of the bulk geometry, encoded in 
the perturbation of E ab . This term can be seen as an external source for perturbations, 
absent in General Relativity. If one regards E ab as an energy-momentum tensor of an 
additional fluid (called the Weyl- fluid), its evolution is connected to the energy density of 
matter on the brane, as one can see from eq. (|42jl. If one neglects the anisotropic stress 
of the Weyl-fluid, then at low energy and superhorizon scales, it decays as radiation, 
i.e. 5E a i, oc a~ 4 . However, the bulk gravitational field exerts an anisotropic stress 
onto the brane, whose time-evolution cannot be obtained by considering the projected 
equations on the brane alone |43j . Rather, the full five-dimensional equations have to be 
solved, together with the junction conditions. The full evolution of E ab in the different 
cosmological eras is currently not understood. However, as we will discuss below, partial 
results have been obtained for the case of a de Sitter brane, which suggest that E ab does 
not change the spectrum of scalar perturbations. It should be noted however, that the 




(45) 



(46) 




(47) 



(48) 
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issue is not settled and that it is also not clear if the subsequent cosmological evolution 
during radiation and matter era leaves an imprint of the bulk gravitational field in the 
anisotropies of the microwave background radiation [IJj • With this in mind, we will, for 
scalar perturbations, neglect the gravitational backreaction described by the projected 
Weyl tensor. 

Considering scalar perturbations for the moment, the perturbed line element on 
the brane has the form 

ds 2 = -(1 + 2A)dt 2 + 2d i Bdtdx i + ((1 - 2ip)S i:j + D ij E)dx i dx j , (49) 

where the functions A, B, E and ip depend on t and x l . 

An elegant way of discussing scalar perturbations is to make use of of the gauge 
invariant quantity 



C = i) + H 5 JL. (50) 

In General Relativity, the evolution equation for ( can be obtained from the energy- 
conservation equation It reads, on large scales, 



H 

C = 1 <$Pnad, (51) 

P + P 

where <5p na d = Sp to t — c 2 5p is the non-adiabatic pressure perturbation. The energy 
conservation equation, however, holds for the Randall-Sundrum model as well. 
Therefore, eq. ()51|) is still valid for the brane world model we consider. For inflation 
driven by a single scalar field 5p na ,d vanishes and therefore ( is constant on superhorizon 
scales during inflation. Its amplitude is given in terms of the fluctuations in the scalar 
field on spatially flat hypersurfaces: 

c = m (5 2 , 



The quantum fluctuation in the (slow-rolling) scalar field obey {(S(j)) 2 ) ~ (H/2ti) 2 , as 
the Klein-Gordon equation is not modified in the brane world model we consider. The 
amplitude of scalar perturbations is [5U| = 4(£ 2 )/25. Using the slow-roll equations 
and eq. (jo^|) one obtains 



A\ 



( 5127T > 


YL 


■20- + V- 


3 




V' 2 


2o 





\k=aH 



(53) 



Again, the corrections are contained in the terms in the square brackets. For a given 
potential the amplitude of scalar perturbations is enhanced compared to the prediction 
of General Relativity. 

The arguments presented so far suggest that, at least for scalar perturbations, 
perturbations in the bulk spacetime are not important during inflation. This, however, 
might not be true for tensor perturbations, as gravitational waves can propagate into 
the bulk. For tensor perturbations, a wave equations for a single variable can be 
derived jSTj. The wave equation can be separated into a four-dimensional and a five- 
dimensonal part, so that the solution has the form hij = A(y)h(x tl )eij, where e%j is a 
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(constant) polarization tensor. One finds that the amplitude for the zero mode of tensor 
perturbation is given by [51] 

A t = ^T7TH 2 F 2 {H/p)\ k=aH , (54) 



25ttM* 



with 



Fix) 



a/1 + x 2 — x 2 sinh 1 [ — 



,x 

where we have defined 



(55) 



TT / O \ 1/2 

7 = (s?) HM «- (56) 

It can be shown that modes with m > 3H/2 are generated but they decay during 
inflation. Thus, one expects in this scenario only the massless modes to survive until 
the end of inflation [SI], [521 ■ 

From eqns. (jB^j) and ([SH]) one sees that the amplitudes of scalar and tensor 
perturbations are enhanced at high energies. However, scalar perturbations are more 
enhanced than tensors. Thus, the relative contribution of tensor perturbations will be 
suppressed, if inflation is driven at high energies. 

Finally, we would like to mention that there are also differences between General 
Relativity and the brane world model we consider for the prediction of two-field brane 
inflation. Usually correlations are separated in adiabatic and isocurvature modes for 
two-field inflation 53J- In the Randall-Sundrum model, this correlation is suppressed 
if inflation is driven at high energies [51]. This implies that isocurvature and adiabatic 
perturbations are uncorrelated, if inflation is driven at energies much larger than the 
brane tension. 

Coming back to cosmological perturbations, the biggest problem is that the 
evaluation of the projected Weyl tensor is only possible for the background cosmology. 
As soon as one tries to analyse the brane cosmological perturbations, one faces the 
possibility that the E 0i terms might not vanish. In particular this means that the 
equation for the density contrast 5 = Sp/p, which is given by {w m = p/p, k is the 
wavenumber) 

k 2 

5 + (2 - 3u m )HS - Qu m {H 2 + H)5 = (1 + co m )5R 00 - uj m —S, (57) 

cr 

cannot be solved as 5i?oo involves 5Eo and can therefore not be deduced solely from 
the brane dynamics 



2.3. Final Remarks on the Randall-Sundrum Scenario 

The Randall-Sundrum model discussed in this section is the simplest brane world 
model. We have not discussed other important conclusions one can draw from the 
modifications of Friedmann's equation, such as the evolution of primordial black holes 
[55], hs connection to the AdS/CFT correspondence (see e.g. [51]- [HI]) and inflation 
driven by the trace anomaly of the conformal field theory living on the brane (see e.g. 
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[H2]-jnS])- These developments are important in many respects, because they give not 
only insights about the early universe but gravity itself. They will not be reviewed here. 

3. Including a Bulk Scalar Field 

In this section we are going to generalize the previous results obtained with an empty 
bulk. To be specific, we will consider the inclusion of a scalar field in the bulk. As we 
will see, one can extend the projective approach wherein one focuses on the dynamics 
of the brane, i.e. one studies the projected Einstein and the Klein-Gordon equation 
jnn], [HZI- As in the Randall-Sundrum setting, the dynamics do not closed, as bulk 
effects do not decouple. We will see that there are now two objects representing the 
bulk back- reaction: the projected Weyl tensor and the loss parameter A$ 2 . I n the 
case of homogeneous and isotropic cosmology on the brane, the projected Weyl tensor 
is determined entirely up to a dark radiation term. Unfortunately, no information on 
the loss parameter is available. This prevents a rigorous treatment of brane cosmology 
in the projective approach. 

Another route amounts to studying the motion of a brane in a bulk space-time. 
This approach is successful in the Randall-Sundrum case thanks to Birkhoff's theorem 
which dictates a unique form for the metric in the bulk [3H] . In the case of a bulk scalar 
field, no such theorem is available. One has to resort to various ansatze for particular 
classes of bulk and brane scalar potentials (see e.g. |SH]-|ZS1)- We will come back to 
this in section 4. 

3.1. BPS Backgrounds 

3.1.1. Properties of BPS Backgrounds As the physics of branes with bulk scalar fields 
is pretty complicated, we will start with a particular example where both the bulk and 
the brane dynamics are fully under control [ZZj (see also [IB] and [EI]). We specify the 
bulk Lagrangian as 



where V(<p) is the bulk potential. The boundary action depends on a brane potential 



where Ub(4>o) is evaluated on the brane. The BPS backgrounds arise as particular case of 
this general setting with a particular relationship between the bulk and brane potentials. 
This relation appears in the study of iV = 2 supergravity with vector multiplets in the 
bulk. The bulk potential is given by 



S = ^| / d 5 x^g~ 5 (J2-| (m 2 + 



(58) 




(59) 




(60) 
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where W(<f>) is the superpotential. The brane potential is simply given by the 
superpotential 

U B = W (61) 

We would like to mention, that the last two relations have been also used in order to 
generate bulk solutions without necessarily imposing supersymmetry [23, (23- Notice 
that the Randall- Sundrum case can be retrieved by putting W = est. Supergravity puts 
further constraints on the superpotential which turns out to be of the exponential type 

E3 

W = Ake^ (62) 

with a = — l/\/l2, l/v3. In the following we will often choose this exponential potential 
with an arbitrary a as an example. The actual value of a does not play any role and 
will be considered generic. 

The bulk equations of motion comprise the Einstein equations and the Klein-Gordon 
equation. In the BPS case and using the following ansatz for the metric 

ds 2 = a{yfr] lxv dx^dx 1 ' + dy 2 , (63) 

these second order differential equations reduce to a system of two first order differential 
equations 

a - = - f . (64) 
a 4 

dW 

Notice that when W = est one easily retrieves the exponential profile of the Randall- 
Sundrum model. 

An interesting property of BPS systems can be deduced from the study of the 
boundary conditions. The Israel junction conditions reduce to 

-\b = ~\b (65) 



and for the scalar field 

At\ 9W \ 

<P\B = -TT-\B 



(66) 



This is the main property of BPS systems: the boundary conditions coincide with 
the bulk equations, i.e. as soon as the bulk equations are solved one can locate the 
BPS branes anywhere in this background, there is no obstruction due to the boundary 
conditions. In particular two-brane systems with two boundary BPS branes admit 
moduli corresponding to massless deformations of the background. They are identified 
with the positions of the branes in the BPS background. We will come back to this later 
in section 5. 

Let us treat the example of the exponential superpotential. The solution for the 
scale factor reads 

a = (1 - Aka 2 x 5 ) 1/4a \ (67) 



Cosmology and Brane Worlds: A Review 



15 



and the scalar field is given by 

0= — Wl -Aka 2 x b ). (68) 
a 

For a — > 0, the bulk scalar field decouples and these expressions reduce to the Randall- 
Sundrum case. Notice a new feature here, namely the existence of singularities in the 
bulk, corresponding to 

a(x 5 )U = (69) 

In order to analyse singularities it is convenient to use conformal coordinates 

du = (70) 
a(x 5 ) 

In these coordinates light follows straight lines u = ±£. It is easy to see that the 
singularities fall in two categories depending on a. For a 2 < 1/4 the singularity is 
at infinity = oo. This singularity is null and absorbs incoming gravitons. For 
a 2 > 1/4 the singularity is at finite distance. It is time-like and not wave- regular, i.e. 
the propagation of wave packets is not uniquely defined in the vicinity of the singularity. 
For all these reasons these naked singularities in the bulk are a major drawback of brane 
models with bulk scalar fields [80]. In the two-brane case the second brane has to sit in 
front of the naked singularity. 

3.1.2. de Sitter and anti de Sitter Branes Let us modify slightly the BPS setting by 
detuning the tension of the BPS brane. This corresponds to adding or substracting 
some tension compared to the BPS case 

U B = TW (71) 

where T is real number. Notice that this modification only affects the boundary 
conditions, the bulk geometry and scalar field are still solutions of the BPS equations 
of motion. In this sort of situation, one can show that the brane does not stay static. 
For the detuned case, the result is either a boosted brane or a rotated brane. We will 
soon generalize these results so we postpone the detailed explanation to later. Defining 
by u(x^) the position of the brane in conformal coordinates, one obtains 

(du) 2 = (72) 

The brane velocity vector d^u is of constant norm. For T > 1, the brane velocity vector 
is time-like and the brane moves at constant speed. For T < 1 the brane velocity vector 
is space-like and the brane is rotated. Going back to a static brane, we see that the bulk 
geometry and scalar field become x^ dependent. In the next section we will find many 
more cases where branes move in a static bulk or equivalently, a static brane borders a 
non-static bulk. 

Let us now conclude this section by studying the brane geometry when T > 1. In 
particular one can study the Friedmann equation for the induced bulk factor 
T 2 - 1 

H 2 = ^7—W 2 , (73) 
16 
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where W is evaluated on the brane. Of course we obtain the fact that cosmological 
solutions are only valid for T > 1. Now in the Randall- Sundrum case W = 4k leading 
to 

H 2 = (T 2 - l)k 2 . (74) 

In the case T > 1 the brane geometry is driven by a positive cosmological constant. This 
is a de Sitter brane. When T < 1 the cosmological constant is negative, corresponding 
to an AdS brane. We are going to generalize these results by considering the projective 
approach to the brane dynamics. 

3.2. Bulk Scalar Fields and the Projective Approach 

3.2.1. The Friedmann Equation We will first follow the traditional coordinate 
dependent path. This will allow us to derive the matter conservation equation, the 
Klein-Gordon and the Friedmann equations on the brane. Then we will concentrate on 
the more geometric formulation where the role of the projected Weyl tensor will become 
transparent |86J, t 87.. Again, in this subsection we will put k 5 = 1. 

We consider a static brane that we choose to put at the origin x 5 = 0. and impose 
6(0, t) = 1. This guarantees that the brane and bulk expansion rates 

AH = d T ^\o, 3H B = ^ T ^/=gl\ (75) 

coincide. We have identified the brane cosmic time dr = ab\odt. We will denote by 
prime the normal derivative d n = -^\od X5 . Moreover we now allow for some matter to 
be present on the brane 

TZ maUer = (- P m,Pm,Pm,Pm). (76) 

The bulk energy-momentum tensor reads 

T ab = ~ (d a( pd b( p) - ^g ab ((d<pf + V) . (77) 

The total matter density and pressure on the brane are given by 

3 3 

P = Pm + ~U B , P = Pm- ~U B - (78) 

The boundary condition for the scalar field is unchanged by the presence of matter on 
the brane. 

The (05) Einstein equation leads to matter conservation 

p m = -3H(p m + p m ). (79) 
By restricting the (55) component of the Einstein equations we obtain 



H 2 = — — —Q — —E + ^— (80) 
36 3^ 9 a 4 



in units of k 2 ,. The last term is the dark radiation, whose origin is similar to the 
Randall-Sundrum case. The quantity Q and E satisfy the differential equations 

Q + AHQ = HTl 
E + AHE = -pT 5 °. 
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These equations can be easily integrated to give 

= & + ^ - — A [dA^ -20)-— A ttraW^, (81) 
36 12 16a 4 J dr yv 1 12a 4 J H dr ' v ; 

up to a dark radiation term and we have identified 

^HfH- (82) 

This is the Friedmann equation for a brane coupled to a bulk scalar field. Notice that 
retarded effects springing from the whole history of the brane and scalar field dynamics 
are present. In the following section we will see that these retarded effects come from the 
projected Weyl tensor. They result from the exchange between the brane and the bulk. 
Notice, that Newton's constant depends on the value of the bulk scalar field evaluated 
on the brane (<f> = <fi(t, y = 0)): 

8nG N (<j>o) = K\U B {<f> Q ) 

3 12 1 ' 

On cosmological scale, time variation of the scalar field induce a time variations of 
Newton's constant. This is highly constrained experimentally (HEMES], leading to tight 
restrictions on the time dependence of the scalar field. 

To get a feeling of the physics involved in the Friedmann equation, it is convenient 
to assume that the scalar field is evolving slowly on the scale of the variation of the scale 
factor. Neglecting the evolution of Newton's constant, the Friedmann equation reduces 
to 

H2= S*G^> pm+ U_t (M) 
3 Hm 8 16 v ; 

Several comments are in order. First of all we have neglected the contribution due to 
the p 2 m term as we are considering energy scales below the brane tension. Then the 
main effect of the scalar field dynamics is to involve the potential energy U and the 
kinetic energy <fi 2 . Although the potential energy appears with a positive sign we find 
that the kinetic energy has a negative sign. For an observer on the brane this looks like 
a violation of unitarity. We will reanalyse this issue in section 5, when considering the 
low energy effective action in four dimensions and we will see that there is no unitarity 
problem in this theory. The minus sign for the kinetic energy is due to the fact that one 
does not work in the Einstein frame where Newton's constant does not vary, a similar 
minus sign appears also in the effective four-dimensional theory when working in the 
brane frame. 

The time dependence of the scalar field is determined by the Klein-Gordon equation. 
The dynamics is completely specified by 

4> + 4H4>+-(-- ou m )p m -^- = —qT + A$ 2 , (85) 

where p m = u m p m . We have identified 

dU B d 2 U B . 

A$ 2 = 0"|o- | . (86) 
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This cannot be set to zero and requires the knowledge of the scalar field in the vicinity 
of the brane. When we discuss cosmological solutions below, we will assume that this 
term is negligible. 

The evolution of the scalar field is driven by two effects. First of all, the scalar field 
couples to the trace of the energy momentum tensor via the gradient of Ub- Secondly, 
the field is driven by the gradient of the potential U, which might not necessarily vanish. 

3.2.2. The Friedmann equation vs the projected Weyl tensor We are now coming back 
to the origin of the non-trivial Friedmann equation. Using the Gauss-Codazzi equation 
one can obtain the Einstein equation on the brane |6b'| . [oT] 

Gab = -^Uh ab + -£-r ab + iTab + ^d a <pd b <p - j^(d<p) 2 h ab - E ab . (87) 

Now the projected Weyl tensor can be determined in the homogeneous and isotropic 
cosmology case. Indeed only the E 00 component is independent. Using the Bianchi 
identity D a G a b = where D a is the brane covariant derivative, one obtains that 

E 00 + AHE 00 = d T (^0 2 + ^Uj + \Hft + ^ Pm (88) 

leading to 

£00 = ^ / dra* (d T (A^ 2 + S -U) + + ^ Pm ) (89) 

Upon using 

Goo = 3H 2 (90) 

one obtains the Friedmann equation. It is remarkable that the retarded effects in the 
Friedmann equation all spring from the projected Weyl tensor. Hence the projected 
Weyl tensor proves to be much richer in the case of a bulk scalar field than in the empty 
bulk case. 

3.2.3. Self-Tuning and Accelerated Cosmology The dynamics of the brane is not closed, 
it is an open system continuously exchanging energy with the bulk. This exchange is 
characterized by the dark radiation term and the loss parameter. Both require a detailed 
knowledge of the bulk dynamics. This is of course beyond the projective approach where 
only quantities on the brane are evaluated. In the following we will assume that the 
dark radiation term is absent and that the loss parameter is negligible. Furthermore, 
we will be interested in the effects of a bulk scalar field for late-time cosmology (i.e. 
well after nucleosynthesis) and not in the case for inflation driven by a bulk scalar field 
(see e.g. |EI|-|S5|)- 

Let us consider the self-tuned scenario as a solution to the cosmological constant 
problem. It corresponds to the BPS superpotential with a = 1. In that case the potential 
U — for any value of the brane tension. The potential U — can be interpreted 
as a vanishing of the brane cosmological constant. The physical interpretation of 
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the vanishing of the cosmological constant is that the brane tension curves the fifth 
dimensional space-time leaving a flat brane intact. Unfortunately, the description of the 
bulk geometry in that case has shown that there was a bulk singularity which needs to 
be hidden by a second brane whose tension is fine-tuned with the first brane tension. 
This reintroduces a fine-tuning in the putative solution to the cosmological constant 
problem [TH]. 

Let us generalize the selftuned case to a ^ 1, i.e. Ub = TW, T > 1 and W is the 
exponential superpotential. The resulting induced metric on the brane is of the FRW 
type with a scale factor 

/ t v l/3+l/6a 2 

a(t)=ao - (91) 



leading to an acceleration parameter 

• = 1^-1 ^ 

For the supergravity value a = — this leads to go = —4/7. This is in coincidental 
agreement with the supernovae results. This model can serve as a brane quintessence 
model [77],jHn]. We will comment on the drawbacks of this model later. See also [HI] 
and [02] for similar ideas. 

3.2-4- The brane cosmological eras Let us now consider the possible cosmological 
scenarios with a bulk scalar field [HEJ^jEZl- We assume that the potential energy of 
the scalar field U is negligible throughout the radiation and matter eras before serving 
as quintessence in the recent past. 

At very high energy above the tension of the brane the non-conventional cosmology 
driven by the (? m term in the Friedmann equation is obtained. Assuming radiation 
domination, the scale factor behaves like 

/ t \ 1/4 

a = a - (93) 



+ 0Mr-) ( 94 ) 



and the scalar field 

.to- 
rn the radiation dominated era, no modification is present, provided 

</> = <j>i (95) 

which is a solution of the Klein-Gordon equation as the trace of the energy-momentum 
of radiation vanishes (together with a decaying solution, which we have neglected). In 
the matter dominated era the scalar field evolves due to the coupling to the trace of 
the energy-momentum tensor. This has two consequences. Firstly, the kinetic energy 
of the scalar field starts contributing in the Friedmann equation. Secondly, the effective 
Newton constant does not remain constant. The cosmological evolution of Newton's 
constant is severely constrained since nucleosynthesis [EE]>t89j. This restricts the possible 
time variation of 6. 
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In order to be more quantitative let us come back to the exponential superpotential 
case with a detuning parameter T. The time dependence of the scalar field and scale 
factor become 



t \ 3 45™ 



(96) 



Mil 

where t e and a e are the time and scale factor at matter-radiation equality. Notice the 
slight discrepancy of the scale factor exponent with the standard model value of 2/3. 
The redshift dependence of the Newton constant is 

G N (z) = / z + 1 \ 4q2 / 5 
G N (z e ) \z e + lj 

For the supergravity model with a = — and z e ~ 10 3 this leads to a decrease by 
(roughly) 37% since nucleosynthesis. This is marginally compatible with experiments 

MM- 

Finally let us analyse the possibility of using the brane potential energy of the scalar 
field U as the source of acceleration now. We have seen that when matter is negligible 
on the brane, one can build brane quintessence models. We now require that this occurs 
only in the recent past. As can be expected, this leads to a fine-tuning problem as 

M 4 ~ Pc (97) 

where M 4 = (T — 1) |^ is the amount of detuned tension on the brane. Of course this is 
nothing but a reformulation of the usual cosmological constant problem. Provided one 
accepts this fine-tuning, as in most quintessence models, the exponential model with 
a = — :^is a cosmological consistent quintessence model with a five dimensional origin. 

3.3. Brief summary 

The main difference between a brane world model with a bulk scalar field and the 
Randall-Sundrum model is that the gravitational constant becomes time-dependent. 
As such it has much in common with scalar-tensor theories j^D] > but there are important 
differences due to the projected Weyl tensor E^ v and its time-evolution. The bulk scalar 
field can play the role of the quintessence field, as discussed above, but it could also 
play a role in an inflationary era in the very early universe (see e.g. |HI]-jHS])- I n any 
case, the cosmology of such a system is much richer and, because of the variation of the 
gravitational constant, more constrained. It remains to be seen if the bulk scalar field 
can leave a trace in the CMB anisotropies and Large Scale Structures (for first results 
see |H7j). 



4. Moving Branes in a Static Bulk 



So far, we were mostly concerned with the evolution of the brane, without referring to 
the bulk itself. In fact, the coordinates introduced in eq. (|14|) are a convenient choice 
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for studying the brane itself, but when it comes to analysing the bulk dynamics and its 
geometry, these coordinates are not the best choice. We have already mentioned the 
extended Birkhoff theorem in section 2. It states that for the case of a vacuum bulk 
spacetime, the bulk is necessarily static, in certain coordinates. A cosmological evolving 
brane is then moving in that spacetime, whereas for an observer confined on the brane 
the motion of the brane will be seen as an expanding (or contracting) universe. In the 
case of a scalar field in the bulk, a similar theorem is unfortunately not available, which 
makes the study of such systems much more complicated. We will now discuss these 
issues in some detail, following in particular [73] and [71] . 

4-1. Motion in AdS-Schwarzschild Bulk 

We have already discussed the static background associated with BPS configurations 
(including the Randall-Sundrum case) in the last section. Here we will focus on other 
backgrounds for which one can integrate the bulk equations of motion. Let us write the 
following ansatz for the metric 

ds 2 = -A 2 (r)dt 2 + B 2 (r)dr 2 + i? 2 (r)d£ 2 (98) 

where dT? is the metric on the 3d symmetric space of curvature q = 0, ±1. In general, 
the function A, B and R depend on the type of scalar field potential. This is to 
be contrasted with the case of a negative bulk cosmological constant where Birkhoff's 
theorem states that the most general solution of the (bulk) Einstein equations is given 
by A 2 = /, B 2 = 1/f and R = r where 

/(O = (99) 



We have denoted by I = 1/k — y— 6/(A 5 k|) the AdS scale and fi the black hole mass 
(see section 2). This solution is the so-called AdS-Schwarzschild solution. 

Let us now study the motion of a brane of tension T/l in such a background. 
The equation of motion is determined by the junction conditions. The method will 
be reviewed later when a scalar field is present in the bulk. The resulting equation of 
motion for a boundary brane with a Z2 symmetry is 

[j 2 + f{r)f /2 = jr (100) 

v 

for a brane located at r [HE]- Here r is the velocity of the brane measured with the 
proper time on the brane. This leads to the following Friedmann equation 

So the brane tension leads to an effective cosmological constant (T 2 — l)// 2 . The 
curvature gives the usual term familiar from standard cosmology while the last term 
is the dark radiation term whose origin springs from the presence of a black-hole in the 
bulk. At late time the dark radiation term is negligible for an expanding universe, we 
retrieve the cosmology of a FRW universe with a non- vanishing cosmological constant. 
The case T = 1 corresponds of course to the Randall-Sundrum case. 
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4-2. Moving branes 

Let us now describe the general formalism, which covers the case of the AdS- 
Schwarzschild spacetime mentioned above. 

Consider a brane embedded in a static background. It is parametrized by the 
coordinates X A (x^) where A = ... 4 and the x^ 1 are world volume coordinates. Locally 
the brane is characterized by the local frame 

4 - ^ (102) 
which are tangent to the brane. The induced metric is given by 

V = gABefe* (103) 
and the extrinsic curvature 

= e A e*D A n B , (104) 
where n A is the unit vector normal to the brane defined by (up to a sign ambiguity) 

g AB n A n B = 1, n A e A = 0. (105) 

For a homogeneous brane embedded in the spacetime described by the metric (JHBJ), we 
have T = T(r), r = r(r) where r is the proper time on the brane. The induced metric 
is 

ds 2 B = -dr 2 + R 2 {r)dY? . (106) 
The local frame becomes 

e A = (T,f, 0,0,0), e A = (0,0,5 A ), (107) 
while the normal vector reads 



n A = (ABr, -By/l + r 2 , 0, 0, 0). (108) 
The components of the extrinsic curvature tensor can found to be 



K i:i = - R&Sij, (109) 

1 d 



K TT = - r51 -(AVl + B 2 f 2 ). (110) 

AB dr 



The junction conditions are given by 

k\( 1 



K, v = —± [t, u - -rh^) . (Ill) 
This implies that the brane dynamics are specified by the equations of motion 



VI + B 2 r 2 R' _ 4 
B ~R ~ "6 

and 



p (112) 



1 d -(AVl + B 2 f 2 ) = -^(2p + 3p), (113) 



AB dr v ' 6 
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where we have assumed a fluid description for the matter on the brane. These two 
equations determine the dynamics of any brane in a static background. 

Let us now close the system of equations by stating the scalar field boundary 
condition 

n A ^ = ^(p-3p), (114) 
where the coupling to the brane is defined by the Lagrangian 

Sbrane = J d 4 xC[tp m , h^], (115) 

where ipm represents the matter fields and 

V = V- ( 116 ) 
This reduces to 

1 -p + 3p). (117) 



2 VI + B 2 r 2 d 

Combining the junction conditions leads to the conservation equation 

p + 3H(p + p) = (p-3p)i (118) 

This is nothing but the conservation of matter in the Jordan frame defined by h^. 

We now turn to a general analysis of the brane motion in a static bulk. To do that 
it is convenient to parametrized the bulk metric slightly differently 

dv 2 

ds 2 = -f 2 (r)h(r)dt 2 + 77- + r 2 dY?. (119) 

h(r) 

Now, the Einstein equations lead to (redefining — > f^-(/> and V — > f^-V^) 



^("+^-«) = -"!(^ + ^) (120) 



3 A . rh' \ 2 IbxjP 



r 2 V 2 * / 
and the Klein Gordon equation 

^ + {v + T +h r = w (122) 

Subtracting eq. ()120|) from (|121|) and solving the resulting differential equation, we 
obtain 

/ = exp J drrA . (123) 

It is convenient to evaluate the spatial trace of the projected Weyl tensor. This is 
obtained by computing both the bulk Weyl tensor and the vector normal to the moving 
brane. With A = y/hf, R = r, B = 1/y/h, this gives 
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This is the analogue of the dark radiation term for a general background. The equations 
of motion can be cast in the form 



K? , kr 2 



ft =-y(^-— W>' 2 ' ( 125 ) 

H' + 4^= -%U-^\ (126) 



K\y =6H+-rH'-3^, (127) 



3 <U> Q ft 

-rn — 3- 

4 r 

where we have defined 

H = ^±. (128) 

This allows to retrieve easily some of the previous solutions. Choosing (ft to be constant 

leads to / = 1, fi is constant and 

1 u 
h — 

P r 4 

This is the AdS-Schwarzschild solution. 

For q = the equations of motion simplify to 

ir<f>> =-^, (130) 



^ = ~T2+^ ( 129 ) 



3 

d(- 

\ft 2 / /' 



-d (^ 4 ) , (131) 



2V (132) 



3 dfx d /H^ 

6 4«;| d<p dcj) \ fi 

In this form it is easy to see that the dynamics of the bulk are completely integrable. 
First of all the solutions depend on an arbitrary function //(</>) which determines the 
dynamics. Notice that 

/ = ^ (133) 
ft 

where /io is an arbitrary constant. The radial coordinate r is obtained by simple 
integration of eq. ()130|) 

a 

r = r e"^^ # . (134) 
Finally the rest of the metric follows from 

h = -^rV fdA4f^ (135) 
3 J dfx 

The potential V then follows (J132)) . This is remarkable and shows why Birkhoff's 
theorem is not valid in the presence of a bulk scalar field. Moreover, it is intriguing 
that the generalization of the dark energy term dictates the bulk dynamics completely. 
It is interesting to recast the Friedmann equation in the form 

H 2 = H + |/ry (136) 
36 
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where H is the Hubble parameter on the brane in cosmic time. One can retrieve 

du 
Id, 



standard cosmology by studying the dynamics in the vicinity of a critical point 7^7 = 0. 



Parametrizing 



/i=^y + 50 2 (137) 
leads to the Friedmann equation 

tf 2 = !(p 2 -^ 2 + 4+°( fl - 4 ) ( l38 ) 

36 Ol 

Here 9 is an arbitrary integration constant. Notice that this is a small deviation from 
the Randall- Sundrum case as 

= r~ B/A (139) 

goes to zero at large distances. Hence, standard cosmology is retrieved at low energy 
and long distance. 



5. Cosmology of a Two— Brane System 

In this section we will once more include an ingredient suggested by particle physics 
theories, in particular M-theory. So far we have assumed that there is only one brane 
in the whole space-time. According to string theory, there should be at least another 
brane in the bulk. Indeed, in heterotic M-theory these branes are the boundaries of the 
bulk spacetime [H|. However, even from a purely phenomenological point of view there 
is a reason to include a second brane: the bulk singularity (or the AdS horizon). As we 
have seen in section 3, the inclusion of a bulk scalar field often implies the presence of a 
naked singularity located away from the positive tension brane. The second brane which 
we include now should shield this singularity, so that the physical spacetime stretches 
between the two branes. Another motivation is the hierarchy problem. Randall and 
Sundrum proposed a two brane model (one with positive and one with negative tension), 
embedded in a five-dimensional AdS spacetime. In their scenario the standard model 
particles would be confined on the negative tension brane. As they have shown, in 
this case gravity is weak due to the warping of the bulk spacetime. However, as will 
become clear from the results in this section, in order for this model to be consistent 
with gravitational experiments, the interbrane distance has to be fixed [01]. This can 
be achieved, for example, with a bulk scalar field. As shown in jM] and j^S], gravity 
in the two-brane model of Randall-Sundrum is described by a scalar-tensor theory, in 
which the interbrane-distance, called radion, plays the role of a scalar field. The bulk 
scalar field will modify the Brans-Dicke parameter (see (HE] and jOZ]) of the scalar field 
and will introduce a second scalar field in the low-energy effective theory, so that the 
resulting theory at low energy in the case of two branes and a bulk scalar field is a 
bi-scalar-tensor theory [98 ,[99j. 

In the following we will investigate the cosmological consequences when the distance 
between the branes is not fixed (for some aspects not covered here see e.g. |lUUj - |107j ). 
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Motivation for this comes, for example, from a recent claim that the fine-structure 
constant might slowly evolve with time [108J. 

5.1. The low-energy effective action 

In order to understand the cosmology of the two-brane system, we derive the low-energy 
effective action by utilizing the moduli space approximation. From the discussion in 
section 3 and section 4 it becomes clear, that the general solution of the bulk Einstein 
equations for a given potential is difficult to find. The moduli space approximation gives 
the low- energy-limit effective action for the two brane system, i.e. for energies much 
smaller than the brane tensions. 

In the static BPS solutions described in the section 3, the brane positions can be 
chosen arbitrarily. In other words, they are moduli fields. It is expected that by putting 
some matter on the branes, these moduli field become time-dependent, or, if the matter 
is inhomogeneously distributed, space-time dependent. Thus, the first approximation 
is to replace the brane-positions with space-time dependent functions. Furthermore, in 
order to allow for the gravitational zero-mode, we will replace the flat spacetime metric 
n^y with g^ u (x a ). We do assume that the evolution of these fields is slow, which means 
that we neglect terms like (d<p) 3 when constructing the low-energy effective action. 

As already mentioned, the moduli space approximation is only a good 
approximation at energies much less than the brane tension. Thus, we do not recover 
the quadratic term in the moduli space approximation. We are interested in the late time 
effects after nucleosynthesis, where the corrections have to be small. 

Replacing T]^ with g^ v {x a ) in ()63j) and collecting all the terms one finds from the 
5D action after an integration over y: 



with a(y) given by (|H7|. The moduli and a represent the location of the two branes. 
Note that the kinetic term of the field (f) has the wrong sign. This is an artifact of the 
frame we use here. As we will see below, it is possible to go to the Einstein frame with a 
simple conformal transformation, in which the sign in front of the kinetic term is correct 
for both fields. 

In the following we will concentrate on the BPS system with exponential 
superpotential from section 3. Let us redefine the fields according to 




4 Kk 



(140) 



with 




(141) 




(142) 
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with (3 = 2c fj l ; and then 

= Q cosh R, a = Q sinh R. (143) 
A conformal transformation g^ u = Q 2 g fJLV leads to the Einstein frame action: 

12a 2 (dQ) 2 



l + 2a 2 Q 2 

(144) 



SeF ~ 2kn 2 (2a 2 + 1) I d ' xV ^ 



2kK 2 5 {2a 2 + 
6 (dR) 2 



K- 



2a 2 + 1 

Note that in this frame both fields have the correct sign in front of the kinetic terms. 
For a —> (i.e. the Randall-Sundrum case) the Q-field decouples. This reflects the 
fact, that the bulk scalar field decouples, and the only scalar degree of freedom is the 
distance between the branes. One can read off the gravitational constant to be 

16nG = 2kK 2 5 (l + 2a 2 ). (145) 

The matter sector of the action can be found easily: if matter lives on the branes, 
it "feels" the induced metric. That is, the action has the form 

S&^S&^gSP) and Sg> = s£)(* 2> (146) 

where g^ denotes the induced metric on each branes. In going to the Einstein frame 
one gets 

SU = SW(* 1: A 2 (Q,R) gflu ) and S$ = S$(V 2 . B 2 {Q, R)g^), (147) 

where matter now couples explicitely to the fields via the functions A and B, which we 
will give below (neglecting derivative interactions). 

The theory derived with the help of the moduli space approximation has the form 
of a multi-scalar-tensoT theory, in which matter on both branes couple differently to the 
moduli fields. We note, that methods different from the moduli-space approximation 
have been used in the literature in order to obtain the low-energy effective action 
or the resulting field equations for a two-brane system (see in particular |109j- 
|113j ). Qualitatively, the features of the resulting theories agree with the moduli-space 
approximation discussed above. 

In the following we will discuss observational constraints imposed on the parameter 
of the theory. 



5.2. Observational constraints 

In order to constrain the theory, it is convenient to write the moduli Lagrangian in the 
form of a non-linear sigma model with kinetic terms 

7^0W, (148) 

where i — 1, 2 labels the moduli <p l = Q and <fi 2 = R. The sigma model couplings are 
here 

12a 2 16 / d \ 
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Notice the potential danger of the a — > limit, the RS model, where the coupling to Q 
becomes very small. In an ordinary Brans-Dicke theory with a single field, this would 
correspond to a vanishing Brans-Dicke parameter which is ruled out experimentally. 
Here we will see that the coupling to matter is such that this is not the case. Indeed we 
can write the action expressing the coupling to ordinary matter on our brane as 

A = a{4>)f- 1/2 {^ or), B = a{a)r 1/2 {<t>, <r), (150) 

where we have neglected the derivative interaction. 
Let us introduce the parameters 

a Q = d Q In A, a R = d R In A. (151) 

We find that (A = 4/(1 + 2a 2 )) 

A = Q-^icoshR)*, (152) 

leading to 

a 2 A 1 Atanhi? 
a Q = 2~q > a R = 4 • ( 153 ) 

Observations constrain the parameter 

$ = j ij a iaj (154) 

to be less than 10" 3 PH. We obt 

ain therefore a bound on 

4 a 2 tanh 2 i? . . 

9 = ; + sr. (155) 

31 + 2a 2 6(1 + 2a 2 ) V ; 

The bound implies that 

a < 1(T 2 , R < 0.2 (156) 

The smallness of a indicates a strongly warped bulk geometry such as an Anti-de Sitter 
spacetime. In the case a = 0, we can easily interpret the bound on R. Indeed in that 
case 

tanhfl=e- fc(ff -* ) , (157) 

i.e. this is nothing but the exponential of the radion field measuring the distance between 
the branes. We obtain that gravity experiments require the branes to be sufficiently far 
apart. When a^O but small, one way of obtaining a small value of R is for the hidden 
brane to become close from the would-be singularity where a(a) = 0. 

We would like to mention that the parameter 9 can be calculated also for matter 
on the negative tension brane. Then, following the same calculations as above, it can be 
seen that the observational constraint for 9 cannot be satisfied. Thus, if the standard 
model particles are confined on the negative tension brane, the moduli have necessarily 
to be stabilized. In the following we will assume that the standard model particles 
are confined on the positive tension brane and study the cosmological evolution of the 
moduli fields. 
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5.3. Cosmological implications 

The discussion in the last subsection raises an important question: the parameter a has 
to be choosen rather small, in order for the theory to be consistent with observations. 
Similarly the field R has to be small too. The field R is dynamical and one would like 
to know if the cosmological evolution drives the field R to small values such that it is 
consistent with the observations today. Otherwise are there natural initial conditions 
for the field R7 In the following we study the cosmological evolution of the system in 
order to answer these questions. 

The field equations for a homogenous and isotropic universe can be obtained from 
the action. The Friedmann equation reads 

" 2 = ^ <* + » + " ; « + ^ + rf^ 2 + ttW* 2 < 158 > 

where we have defined Q = exp 0. The field equations for R and read 



14- 2rv 2 
R + 3HR= - 8tiG \ 



dV eS dW eS 



+ 3#0= -8ttG 



dR OR 

+ 4 1} ( Pl - 3pi) + (p 2 - 3p 2 )] (159) 
l + 2a 2 



12a 



2 



OVcs | dW eS 



+ 4\ Pl - 3 Pl ) + af\p 2 - 3p 2 )] . (160) 
The coupling parameter are given by 

4 1} =-7^ 4 2) = -7^ (161) 



(1) _ tanhfl (2) _ (tanh_g) 1 

" TT2^' fl ~ l + 2a 2 • [ b } 

We have included matter on both branes as well as potentials V c fj and W e g on each 
branes. We now concentrate on the case where matter is only on our brane. In the 
radiation dominated epoch the trace of the energy-momentum tensor vanishes, so that 
Q and quickly become constant. The scale factor scales like a(t) oc t 1 / 2 . 
In the matter-dominated era, the solution to these equations is given by 

/ a \ -3-2a 2 /3 / I \ 2/3-4a 2 /27 

Pi = Pe( — ) , a = a e ( T ) (163) 

together with 



3 CL e \t e J Vt e , 

as soon as t ^> t e . Note that R indeed decays. This implies that small values of R 
compatible with gravitational experiments are favoured by the cosmological evolution. 
Note, however, that the size of R in the early universe is constrained by nucleosynthesis 
as well as by the CMB anisotropies. A large discrepancy between the values of R 
during nucleosynthesis and now induces a variation of the particle masses, or equivalently 



a e J \t, 



+ -\n-,R = Ro -) +Ri[-) , (164) 
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Newton's constant, which is excluded experimentally. One can show that by putting 
matter on the negative tension brane as well, the field R evolves even faster to zero [98J . 
This behaviour is reminiscent of the attractor solution in scalar-tensor theories j!15j . 

In the five-dimensional picture the fact that R is driven to small values means that 
the negative tension brane is driven towards the bulk singularity. In fact, solving the 
equations numerically for more general cases suggest that R can even by negative, which 
is, in the five-dimensional description meaningless, as the negative tension brane would 
move through the bulk singularity. Thus, in order to make any further progress, one has 
to understand the bulk singularity better|. Of course, one could simply assume that 
the negative tension brane is destroyed when it hits the singularity. A more interesting 
alternative would be if the brane is repelled instead. It was speculated that this could 
be described by some effective potential in the low-energy effective action |98j . 

6. Epilogue: Brane Collision 

We have seen that brane world models are plagued with a singularity problem: the 
negative tension brane might hit a bulk singularity. In that case our description of the 
physics on the brane requires techniques beyond the field theory approach that we have 
followed in this review. It is only within a unified theory encompassing general relativity 
and quantum mechanics that such questions might be addressed. String theory may be 
such a theory. The problem of the nature of the resolution of cosmological singularities in 
string theory is still a vastly unchartered territory. There is a second kind of singularity 
which arises when two branes collide. In such a case there is also a singularity in the 
low energy effective action as one of the extra dimensions shrink to zero size. It was 
speculated that brane collisions play an important role in cosmology, especially in order 
to understand the big bang itself |llfi| - [T2*l] . 

In heteroric M-theory the regime where the distance between the branes becomes 
small corresponds to the regime where the string coupling constant becomes small and 
therefore a perturbative heterotic treatment may be available. In particular for adiabatic 
processes the resulting small instanton transition has been thoroughly studied (see e.g. 
|125j and [126J). Here we would like to present an analysis of such a collision and of 
the possible outcome of such a collision. A natural and intuitive phenomenon which 
may occur during a collision is the existence of a cosmological bounce. Such objects are 
not avalaible in 4d under mild assumptions, and therefore can be exhibited as a purely 
extra dimensional signature. 

We will describe a ci-dimensional theory with a scalar field and gravity whose 
solutions present a cosmological singularity at t = 0. It turns out that this model is 
the low energy approximation of a purely (d + 1) dimensional model where the extra 
dimension is an interval with two boundary branes. The singularity corresponds to 
the brane collision. In the (d + 1) dimensional picture, one can extend the motion 

| For a — the theory is equivalent to the Randall-Svmdrum model. In this case the bulk singularity 
is shifted towards the Anti-de Sitter boundary. 
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of the branes past each other, hence providing a continuation of the brane motion 
after the collision. The (d + 1) dimensional space-time is equivalent to an orbifold 
where the identification between space-time points is provided by a Lorentz boost. 
These spaces are the simplest possible space-times with a singularity. As with ordinary 
spatial orbifolds, one may try to define string theory in such backgrounds and analyse 
the stringy resolution of the singularity. Unfortunately, these orbifolds are not stable in 
general relativity ruling them out as candidate stringy backgrounds. Let us now briefly 
outline some of the arguments. 

We have already investigated the moduli space approximation for models with a 
bulk scalar field. Here we will consider that at low energy the moduli space consists of 
a single scalar field coupled to gravity 

S = J d d xV=g(R-\(d<pfy (165) 
Cosmological solutions with 

ds 2 = a 2 {t)[-dt 2 + dx^Xi] (166) 
can be easily obtained 



a = a(l)\t\&, = 0(l) + 6y 2( ^_ 2 1) In \t\, (167) 

where e = ±1. There are two branches corresponding to t < and t > connected by 
a singularity at t — 0. 

So what is the extra dimensional origin of a such a model? One can uplift the 
previous system to (d + 1) dimensions by defining 

V> = e 7 * (168) 



h» = r 4/{d - 2) 9,„ (169) 



and 



where 7 = y (d — 2)/8(d — 1). Consider now the purely gravitational (d+1) dimensional 
theory with the metric 

ds 2 d+l = ^dw 2 + g^dxV, (170) 

where w G [0,1]. The two boundaries at w = and w = 1 are boundary branes. 
The dimensional reduction on the interval w G [0,1], i.e. integrating over the extra 
dimension, yields the effective action ()165j) provided one restricts the two fields ij){x^) 
and g^x^) to be dependent on ^-dimensions only. 

Let us now consider the nature of -dimensional space-time obtained from 

the solutions ()167j) . The (d + 1) dimensional metric becomes 

ds 2 d+1 = B 2 t 2 dw 2 + i]^dx^dx v (171) 

for a given B depending on the integrations constants 0(1) and a(l). The geometry of 
space-time is remarkably simple. It is a direct product R d ~ l x M where M is the two 
dimensional compactified Milne space whose metric is 

ds 2 M = -dt 2 + BH 2 dw 2 . (172) 
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Using the light cone coordinates 

x ± = ±te ±B2w (173) 
the metric of Milne space reads 

ds 2 M = dx + dx-, (174) 

coinciding with the two dimensional Minkowski metric . There is one subtlety here, the 
original identification of the extra-dimensional interval is here transcribed in the fact 
that Milne space is modded out by the boost 

x ± ^ e ±2B* x ±^ (1?5) 

as we have identified the interval with S l /Z 2 and the boundary branes are the fixed 
points of the Z 2 action as in the Randall- Sundrum model. 

The two boundary branes collide at x = 0, their trajectories are given by 

x% = ±t,xf = ±te ±B \ (176) 

At the singularity one can hope that the branes go past each other and evolve henceforth. 
Unfortunatly, Horowitz and Polchinski have shown that the structure of the orbifold 
space-time is unstable |127j . By considering a particle in this geometry, they showed 
that space-time collapses to a space-like singularity. Indeed one can focus on a particle 
and its n— th image under the orbifold action. In terms of collision the impact parameter 
b becomes constant as n grows while the centre of mass energy y/s grows like coshn5 2 . 
As soon as n is large enough, 

G^fs>b d ~ 2 (177) 

the two particle approach each other within their Schwarzschild radii therefore forming 
a black hole through gravitational collapse. So the orbifold space-time does not make 
sense in general relativity, i.e. not defining a time- dependent background for string 
theory. 

Hence, it seems that the most simple example of brane collision needs to be 
modified in order to provide a working example of singularity with a meaningful string 
theoretic resolution. It would be extremely relevant if one could find examples of stable 
backgrounds of string theory where a cosmological singularity can be resolved using 
stringy arguments. A particularly promising avenue is provided by S-branes where a 
cosmological singularity is shielded by a horizon (see e.g. j2H] and references therein). 
Time will certainly tell which of these approaches could lead to a proper understanding 
of cosmological singularities and their resolutions, an issue highly relevant to brane 
cosmology both in the early universe and the recent past. 



7. Open Questions 

In this article we have reviewed different aspects of brane cosmology in a hopefully 
pedagogical manner reflecting our own biased point of view. Let us finally summarize 
some of the open questions. 



Cosmology and Brane Worlds: A Review 



33 



• In the case of the single brane model by Randall & Sundrum, the homogeneous 
cosmological evolution is well understood. An unsolved issue in this model, however, 
is a complete understanding of the evolution of cosmological perturbations. The 
effects of the bulk gravitational field, encoded in the projected Weyl-tensor, on 
CMB physics and Large Scale Structures are not known. The problem is twofold: 
first, the bulk equation are partial, non-linear differential equations and second, 
boundary conditions on the brane have to be imposed. The current formalism have 
not yet been used in order to tackle these problems (for perturbations in brane 
world theories, see |128j - [Tlo] ; a short review on brane world perturbations is given 
in UHI). 

• Models with bulk scalar fields: Although we have presented some results on the 
cosmological evolution of a homogeneous brane, we assumed that the bulk scalar 
field does not strongly vary around the brane. Clearly, this needs to be investigated 
in some detail through a detailed investigation of the bulk equations, presumably 
with the help of numerical methods. Furthermore, for models with two branes, 
the cosmology has to be explored also in the high energy regime, in which the 
moduli-space approximation is not valid. Some exact cosmological solutions have 
been found in (201 and j!48j . 

• Both the bulk scalar field as well as the interbrane distance in two brane models 
could play an important role at least during some part of the cosmological history. 
Maybe one of the fields plays the role of dark energy. In that case, it is only natural 
that masses of particles vary, as well as other parameter, such as the fine structure 
constant a em |149j . Details of this interesting proposal have yet to be worked out. 

• The bulk singularity, which was thought to be shielded away with the help of 
a second brane, seems to play an important role in a cosmological setting. We 
have seen in section 5 that the negative tension brane moves towards the bulk 
singularity and eventually hits it. Therefore, cosmology forces us to think about 
this singularity, even if it was shielded with a second brane. Cosmologically, the 
brane might be repelled, which might be described by a potential. Alternatively, 
one might take quantum corrections into account in form of a Gauss-Bonnett term 
in the bulk [TSnj-fTKH]. 

• Brane collisions provide a different singularity problem in brane cosmology. String 
theory has to make progress in order to understand this singularity as well. From 
the cosmological point of view, the question is, if a transition between the brane 
collision can provide a new cosmological era and how cosmological perturbations 
evolve before and after the bounce jl "> Ij- l(i() . 

These aspects of brane cosmology raise very interesting questions. It is clear that 
cosmology will continue to play an important role for testing our ideas beyond the 
standard model of particle physics. 
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